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I:  Introduction 


Imagine  you  hove  a  greenhouse  in  Belgium.  Hero  you 
have  a  natural  solar  collector  that  doesn't  cost  you  anythin--', 
out  you  still  spend  thousands  of  dollar:;  a  year  on  heating 
fuel.  What's  core,  and  this  is  the  really  frustrating  part, 
on  the  few  beautiful  days  there  are,  you  have  to  open  win¬ 
dows  to  let  out  the  hot  air  d urine  the  day,  and  then  heat 
with  oil  during  the  night  .  You  think  this  is  crazy,  which 
it  is,  and  you  decide  to  capture  scr.e  of  that  heat  you're 
letting  out  during  the  day,  and  use  it  at  night. 

So  what  do  you  do?  There  are  several  possibilities. 

You  can  put  some  rocks  under  your  greenhouse,  and  then  blow 
the  hot  air  through  durinr  the  day  and  blow  cold  air  through 
the  hot  rocks  during  the  night.  Or  you  can  fill  old  coke 
bottles  with  paraffin,  or  sor.e  other  phase  change  material, 
and  run  the  hot  air  through  them.  That  way  you  have  addi¬ 
tional  heat  capacity  in  the  same  volume,  due  to  the  heat  of 
fusion  of  the  material.  Also  the  storage  will  be  at  a  more 
or  less  constant  temperature,  the  temperature  of  fusion,  while 
it  is  being  loaded  and  unloaded  with  heat.  Another  thing  you 
might  do  would  be  to  put  some  plastic  tubes  under  the  greenhouse, 
put  the  phase  change  material  in  between  the  tubes,  and  then 
blow  air  through  the  tubes.  The  hot  or  cold  air  will  then 
give  or  take  up  heat  as  it  nasses.  Finally,  you  could  put 
a  phase  change  material  in  tubes,  and  blow  the  air  through 
the  tubes.  There  are  other  possibilities,  and  the  different 
geometries  and  encapsulation  methods  are  limited  only  by 
one's  imaginiati on. 

What  would  be  good  to  know  is  how  to  make  a  heat  storage 
system  that,  over  the  lifetime  of  the  greenhouse,  will 
maximize  the  amount  of  money  saved  by  acting,  instead  of 
continuing  to  throw  out  the  heat  captured  during  the  day. 

For  the  storage  system  you'd  like  to  know  how  to  get 
the  most  heat  in  and  out  in  the  shortest  amount  of  time; 


that  will  corrennonii  to  the  nininun  amount  of 


electric  energy 

needed  for  blowing  the  air.  In  fact  for  a  storage  system,  one 
can  name  a  criterion,  a  coefficient  of  performance ,  just  like 
a  heat  punn  has,  which  tell"  how  much,  heat  is  saved  fcr  a  reiver 
quant  it;/-  of  work,  or  punninr  energy. 

The  goal  of  this  report  in  to  use  the  results  of  ur.  ex¬ 
perimental  investigation  three  heat  exchanger  storage  units, 
in  order  to  decide  which  is  best.  the  advantages  and  disad¬ 
vantages  of  the  encapsulation  methods  and  geometries  will  be 
discussed;  together  with  this,  a  theoretical  model  and  its 
relation  to  the  experimentally  observed  results  will  be  pre¬ 
sented  . 

II:  Review  of  Literature 


For  his  doctoral  thesis  at  the  von  Karman  Institute, 

Theun i s sen ( 8 )  has  done  a  review  of  the  state  of  the  art  in 
heat  exchanger  modeling  and  the  chemical  processes  which  are 
involved  in  heat -o f - f u s io n  storage  systems.  As  the  present 
work  is  complementary  to  the  T'neunissen  thesis,  a  short  re¬ 
view  of  some  additional  articles  is  presented  here. 

The  world  expert  for  phase  change  materials  is  Dr. 

Maria  Telkes.  She  first  suggested  the  idea  of  using  Glauber's 
salts,  paraffin,  and  eutectic  salts  for  storing  heat  in  the 
1950's  .  In  Solar  Engineering  (l)  she  notes  in  a  review  of 
the  costs  and  benefits  of  various  storage  systems,  that  the 
people  who  work  with  phase  change  materials  should  be  physical 
chemists,  since  many  of  the  problems  associated  with  these 
materials  are  chemical  problems.  Examples  of  this  are  the 
separation  of  the  liquid  from  the  solid  in  hydrar.ed  salts, 
and  supercooling  during  the  r e so  1 i d i f i ca t ion  of  the  melted 
material.  She  also  points  out  the  need  for  attention  to  be 
paid  to  the  method  of  encapsulation  of  the  material. 

In  the  Journal  of  Solar  Eng i n eer ing ( 2  )  a  numerical  sim¬ 
ulation  is  presented  which  is  similar  to  the  one  done  in  the 
present  program.  That  is,  the  model  consists  of  energy  balances 


k 
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both  on  the  heat  transfer  fluid  and  on  the  solid  ohase, 
which  receives  the  heat.  In  the  trade  publication,  Tolar 
Ape  (  3 ,  *» )  passive  latent  heat  storage  systems  arc  discussed 
as  a  r  e  h  y  b  r  i  d  s  y  3 1  e  r.  s  and  the  need  for  low -  energy -  co  n  s  u  n  i  n  g 
pumoi n g . 

in  another  article  by  ?elkes( 5 ) »  the  need  to  consider 
systems  which  use  electric  resistance  heating  during  off- 
peak  utility  periods,  and  the  need  to  consider  the  encapsulation 
are  discussed.  friith  et  al.  {(1)  note  that  the  limiting 
factor  in  heat  exchanger  nerfonance  of  storage  systems  is 
the  resistance  to  heat  transfer  in  the  solid  nhase.  In 
modeling  the  r.icroscoric  phenomena  of  the  phase  change  heat 
exchange,  they  stress  the  importance  of  natural  convection  in 
determining  the  shame  of  the  meltinr  front  and  in  thus  i  n  - 
fluencinr  the  rate  of  heat  transfer. 

The  heat  transfer  coefficient,  h  ,  models  all  the  heat 
transfer  from  the  fluid  to  the  solid  phase,  or  vice  versa. 

This  is  always  based  on  empirical  correlations,  and  thus  may 
be  a  dominant  cause  of  error  in  any  modeling  of  heat  exchangers. 

In  this  brief  sample  of  the  literature  on  the  latent -heat 
storage  subject,  there  are  three  problems  which  are  continually 
mentioned.  First,  there  is  the  need  to  model  mathematically 
the  performance  of  the  heat  storage  unit.  Fecond,  there  is 
the  need  to  devise  low-cost  en c annul a t i on  methods.  Finally 
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III:  T  hpor y 

We  have  the  first  law  of  t  he  rrnodyn  am  i  ~  s ,  which  says 
that,  for  a  riven  device. 

Hate  of  heat  Hate  of  h r- a t  Rate  of  heat 

going  in  +  Production  -  goinr  out 

For  the  air  going  into  a  heat  accumulator,  one  can 
consider  the  rate  of  heat  going  in  to  be: 

U.  .A _C  .  T.  Convection  of  hot  air  in 

x n t  3  p  ,  a  i  r  in 


The  rate  of  heat  going  out  equals 

P__.  U.  AC  .  T  . 

\air  mt  3  p,air  out 


Convection  of  hot  air  out 


+  hiAi(*air^y^“  TFC’!(y))  Heat  s!ven  to  solid 

from  air 

+  h _ A ,  ( T  .  (v)-  T  ,  .  .  (y))  Heat  given  to  exterior 

2  4  air  ■  ambient  throu,.h  valls  of  storage 


Where  A^  Heat  exchange  surface  area 

A^=  Accumulator  cross-sectional  area 

A^=  Accumulator  interstitial  area 

A^  =  External  surface  area  of  acunulator 

The  rate  of  heat  oroducticn,  in  the  case  of  a  machine 
degrading  electrical  or  mechanical  energy  equals: 


p  V  .  C  .  dT  .  /dt 
V  <»*►*  air  p,air  air 

If  we  put  all  these  terms  into  an  equation  according  to 
the  first  law  of  thermodynamics,  we  have: 

P  .  U.  ,A,C  .  T.  +  P  .  V  .  C  .  dT  .  /dt= 

I  air  int  3  p,air  in  •  air  air  p,air  air 


0  .  U.  AC  .  T  +  h  A  ( T  .  (y)-T  (y)  )  +  h  A ,  (T  .  (y)-T  ,.  (>')  ) 

\  air  int  3  p,air  out  1  1  air  PC..  2  U  air  ambient 


Nov  if  we  take  this  equation  and  divide  all  terms  by 
the  volume  of  the  accumulator,  we  have  (  remember  i  nr.  that 
volume  is  equal  to  length  times  cross-sectional  area): 


Q  .  U.  A_C  .  T.  n  .  Ik  ,A,C  .7  „  Q 
'  air  int  3  p,air  m  pair  int  3  r,air  out  p 


.  V  .  C  .  dT  . 
air  air  n , a  1 r  air 


a2l 


A„L 


h ,  A , ( T_ . J  y ) -T ( y  )  )  V  U  (  ^  i  r (  >*  }  "T~  -  -  -  *  (  -v  >  > 


‘  1  1 4  *  a  i  r 


PCM 


ambient 


a2l 


A„L 

d. 


From  this  equation,  we  let 


6  s  =  A3/A2 
£v  =  Vair/A2L 
as=  Ai/A2L 

heat  transfer  coefficient  between  air  and  solid  inside 


the  heat  accumulator 

heat  transfer  coeffi< 
of  the  heat  accumulator 


h2=  heat  transfer  coefficient  between  air  inside  and  outside 


Rewriting  the  energy  balance  with  the  new  notation, 
we  have: 

p  .  U.  .f  c  .  (t. 

\  air  mt^s  p,air\  m,: 


.air 


1  o  u  t  ,  a  i  r)  ,C  0  C  ■  d  T  •  /  d  t 
_  +tvy'a  p,air  air  = 


M“8(T«ir(-v)-7PCM(=’>>  *  h2Sc(Tair(!’)-Ta»bient(y,> 

Letting  the  length  of  our  accumulator  become  infinitesimal, 
we  have: 

-0  .  U.  C  .  dT  .  /dy  +£  P  .  C  .  d'f  .  /dt  = 

\air  mt^s  p,air  air  Vv\ air  p,air  air 


or  a  complete  statement  for  the  energy  balance  on  the  air  . 


L 
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Similarly,  we  have  an  energy  balance  for  the  PCM,  that 
is  that  all  the  heat  it  receives  results  from  heat  given  uo 
by  the  air: 


PCMVPC? 


'p,PCM  PCM 


/d  t 


liow  we  begin  to  make  ass  unctions.  First  of  all,  let  us 
say  that  the  external  losses  and  the  accumulation  of  heat  in 
the  air  are  negligible  when  c  on par  ed  to  the  convection  terms 
in  the  energy  balance.  The  energy  balance  then  becomes: 


fairUint  sCp, a ir dTair /dy  =  Vs(Tair (  y  )_TPCM(y  }  } 


together  with  the  equation  above  it  on  this  page. 


The  next  assumption  is  that  the  T^,,  is  constant  with 
time  and  space,  at  the  temperature  of  fusion  of  the  PCM, 
while  the  fusion  is  occuring. 

Thus  we  obtain: 


dT  .  (y) 

air  ‘ 


T  .  ( y  )-T 

air  fusion 


h ,  a 
1  s 


0  .  u.  .£  c 

r  air  mt^s  p,air 


dy 


Remembering  the  definitions  of  a^  and^,  we 


have : 


dTair<*> 


T  .  ( v )-T 
air  fusion 


Vi 


dy 


A  •  I/O  .  U  .  C 
3  »air  in t  p,air 


and  we  see  that  the  multiplier  of  dy  equals 


Vi 


m  C  .  L 
p,air 


where  m  is  the  air  mass  flow-rate. 


i  4  - _ 


»ov  i  n  t  egr  at  in  r  both  sides  of  the  energy  balance,  wehave  : 

?  .  (v) 

air*,.,  ,  , 

d  ;  .  ( y  )  r 

air  h . A .  i  v 


. ! •  (y)  -  T  .  ) 

air  fur.io  r. 


mC  .  L  J  C. 
p  ,  a  i  r  'J 


vr.ich  ,  evaluate:,  is: 


In  ( •  (  y  ) 


: us i on 


-  Inf  T.  ~7  .  )  = 

in  t u  s i o  n 


p  ,  a  i  r 


from  which  we  arrive  at: 


T  [  y  J  _  n 
a  i  r  w  ‘fusion 


-  a' 

l n  fusion 


=  expf -St'  y/L) 


where  St  is  defined  as  h  Aj 


p ,  a  lr 

Physically  this  dimensionless  group  the  Stantor.  number 
gives  us  the  ratio  of  the  heat  transfered  in  the  heat  accumulato 
to  the  heat  convected  through  the  accumulator. 

From  this  derivation  we  can  construct  the  following  r  r n  r  h : 


(%)  -  T, 


g  tor i nr 
eduction  a 
o o p  solar 


heat  stor 


Q\.0\I ■>614 


Oi.irrnn  of  Collector- ftorare  Loon 


One  nay  have  noticed  by  now  that  not  much  bar  Veen  sail 
about  h^,  the  heat  transfer  coefficient.  This  :s  a  pro:  >" 
which  can  be  addressed  only  by  means  of  empirical  correlatio 
Ve  knov  that  original  ly,  the  f  tant.cn  number  locally  enualn 

h  1  ^  a  ir  Cv> ,  a  ir  J  i  tit 

NV*  X  Cr,air  X  '  UintDn  =  Su/FrH'' 


To  it  rer.air.3  to  find  the  Nusselt  number  by  means  of  e^riric 
correlations.  For  a  packed  bed,  ve  have: 


*>  ('  7 

nr  rv7  Op,Jpr*  '  ‘ 

Nu*  (  1  +  1  ,5(  »-£  )  )(2+.6f,h  1  +  (  — - “ - “To - - - 

v  1+2. UU( Ft" /a-1 )He 


)  Fr 


While  for  the  tube  cross  flow. 


Nu*  .3  He 

and  for  the  tube  parallel  flow,  in  the  turbulent  rcrirt 

r  r 

Ku*  0.2?  I’r  *  He 


by  ri.'arranciri!'  the  terns  in  the  energy  balance  equations 
:e  arrive  at  the  follow i nr: 


(IT  .  / d  (  y  /  L  ) 

air  '  * 


=  r. 


i A  i !(°~  i 


0  .  !’  .  A 

r  p  ,  a  i  r  :  r.  t  ? 


<17 . /  d  r  =  h. 


£  „  Pair%,,ir 


/0  .  c 

v  a  i  r  ’ 


1  -£ 


P  w-:c 


PC”,  P.PC’ 


These  ?  equations  and  their  discretization  fore  the  subject  of 
the  work  of  Ouestois  (Q). 

A  brief  summary  of  the  numerical  approach  follows.  We  have 
a  backward  time,  backward  space  discretization  scheme: 


dT/dr  =  T(  i  ,j  )-T(  i-1  ,.j  )/ aT  for  i  =  1  to  n 


d  7  /  d ( y / L  )  =  T  (  i  ,  j  )  -  T( i , j - 1 ) /  A  ( y / L )  for  j=  2  to  m 


This  is  straightforward;  the  only  part  which  requires  a  bit 
of  thought  is  what  to  introduce  as  the  heat  capacity  of  the 
PCI!  when  the  fusion  is  occuring.  The  approach  chosen  is  fir; 
to  use  a  delta  function,  of  finite  width  to  express  a  "heat 
capacity  of  fusion".  Normally  this  would  be,  with  a  precise 
fusion  temperature; 


This  is  first  modelled  by  havinp  a  finite-width  spike  centered 
on  the  fusion  temperature; 


T 


and  then  by  a  smoother  distribution  of  heat  capacity  about  t.h< 

a  s 


temperature  of  fusion. 


shown  on  the  ton  of  the  next  ->a  -  <* . 


C,  ?  H 


The  simulation  therefore  incrni:>‘'r.  the  heat  jarici  *.y  of 
the  PC"  as  ti:e  fusion  temperature  is  unnroached,  and  the 
phenomenon  of  the  tor.reraturo  remaining  constant  during 
fusion  is  simulated. 


IV :  Comnar i s on  of  Different  ftornme  Systems 

To  compare  the  different  methods  of  heat  storage,  it  is 
important  to  consider:  a  water  storage,  a  rock  bed  storare, 
a  paraffin  phase  chanre  storage,  and  a  CaCl^H^C1  system: 


all  these 

systems  ar 

e  compared  for 

fluid  . 

Water 

Rock  bed 

1  . 

Pipes 

1 .  Rocks 

First 

2  . 

Storage 

2.  Air  blower 

Cost  s 

tank 

3.  Air  ducts 

3. 

Heat 

4 .  Walls  to 

exchanger 

hold  rocks 

U. 

Blower 

5. 

Water 

pump 

6. 

Ducting 

for  air 

air  as  the  solar- heated 
Paraffin  CaCl  fdigO 


Paraffin 

1  . 

CaCl?6P, 

Encarsu- 

2  . 

Fncnpsul 

Nation 

lati on 

3  : 

Ducting 

Ducting 

Blower 

u . 

Blower 

Volume 

required 


1 .  Water 
Volume 

2 .  Heat 
Exchan  t-;er 
volume 


Volume  of  volume  of  PC'*  +  Encapsulation 

rock  bed 


Oper at i on 
costs 


In  all  cases,  a 
the  pressure  lo 
rate  of  the  air 


function  of  the  heat  transfer  coefficient 
is  across  the  storage,  and  the  volume  flow 


across 


r 


Fror.  this  qual  i  tat :  vc  evaluation  v*»  r,o  *o  a  quantitative 
comparison  of  volume  requ  i  r  aments  for  d  i  f Keren t  storage  r.ed  i  a 
q  i  v  e  r.  a  set  storage  c  a  n  a  c  i  t  y  of  SCO  ncrajoulcs,  »i'/on  a  tor.n- 


eraturo  c  h :in/re 

i'  r  o  - 

16CG  to  ? 5 ° C . 

For  water. 

C  = 

p 

1  Tea  1  /K0-°C 

=  ‘i.io  v;  j  o  u  1  e  /  y.  q 

_Oq 

Thus  the  volume 

o  f  V 

’,*1  t  Af* 

5  X  1  0  8  /  (  .  1  o  x 

*3 

1  0  ' 

X  2  X  1 0  1  )  = 

5.07  X  10U  KG  = 

5^.7  T3 

For  rocks  , 

C  s 
n 

.9  Kjoulc/KG- 

°C,  =  2500  KG/' 

,.3 

ana  we 

will  assume  a  void  fraction  of  C.h.  One  cubic  meter  of 

3 

storage  will  then  contain  0.6  V.  of  rock,  or  1SC0  KG. 

V<e  then  have: 

R  X 

5  X  10/(900  X  20)=  2.3  X  10  KG 

2.8  X  1 0^/(1. 5  X  1  03  )  =  18.5  ",3  of  rock  bed 


For  Paraffin,  C  =  2000  J0ULK/K0-°C  anc  we  see  that 

p 

5  X  10^/(2000  X  20)  =  1.25  X  10 ’  KG  are  required.  Therefore 
the  mass  of  paraffin  required  will  be  lower  than  that  of  water 
or  rocks,  but  one  must  remember  that  encapsulation  is  required, 
which  will  have  void  space  as  well  as  volume  itself. 

For  CaClgfiH^O,  the  heat,  sensible,  of  solid  and  liquid 
phases  is  2000  Joule/KG-°C,  the  heat  of  fusion  is  170000  joule/ K 
the  temperature  of  fusion  is  29°C  ,  the  solid  density  is 
1.5  KG/Liter,  and  the  liquid  density  is  1  .  1*  KG/Liter.  We  have: 

5  X  108/  ( 20°C  X  2000  +  170000)  =  2.U  X  103  KG 

2 . U  X  1 03  /  1500  =  1.6  M3 

As  with  paraffin,  there  is  encapsulation  and  void  space  required 
but  there  is  still  a  significant  volume  reduction  over  rocks 

or  water  . 

Together  with  less  volume,  there  is  a  smaller  bed  through 
which  to  pumn  the  hot  air  from  the  solar  collector,  and  this 
will  result  in  lower  blowing  or  rumrinr  costs, 
less  pressure  loss  with  the  smaller  accumulator 


since  there  is 


V.  K xne r  i  r.p.nt  u  1  Pyston 


As  mentioned  before,  three  heat  accumulators,  or  three 
methods  of  encar su 1  at i on  ,  were  tested.  The  concepts  are: 
tube  cross  flow,  tube  parallel  flow,  and  a  racked  bed.  The; 
are  diagrammed  below. 
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Tube  Parallel  Flow 


Tube  Cross  Flow 


Packed  Bed 


In  addition  to  thermal  characteristics,  such  as  the  Stanton 
number,  there  are  several  other  means  by  which  the  accumulators 


may  be  compared,  as  tabulated  belo-. 


Heat  exchange 
Surface  area 


Tube 

Cross  PI  o’. 

16.1  ”2 


Packed  Bed 

1 1* .  3  M? 


Tube 

Parallel  Flow 

5  •  9  M2 


Void  fraction 


0.58 


0.35 


Vo  1 .  PCM 
Vol.  total 


0.26 


Vol.  e  r.  c  a  r  r.  u  1  a  t  i  o  n 
Vol.  total. 

Heat  Capacity/M^ 

(  1  5°C  to  35 °C ) 


77  m  e  p,  a  j  o  u  1  e 


0 . 03 


1  5  4  mepajoule  1(38  me pa joule 


Total  Volume 


333  Liter: 


175  Liters  210  Liter: 


8  mm  H^O 


1  .  5  nn  K,, 0  3  mm  Ho0 


Pressure  drop 


From  a  first  examination  of  these  data  for  the  three 


accumulators,  one  can  see  some  obvious  tradeoffs.  For  exam¬ 
ple  the  tube  parallel  flov,  while  having  the  highest  heat 
capacity  per  unit  volume,  also  has  the  least  area  for  heat 
exchance.  The  volume  of  encapsulation  is  also  low,  which 
results  in  a  low  first  cost,  but  this  advantage  is  offset  by 
the  low  heat  exchance  surface  area,  which  results  in  higher 
operating  costs  due  t.  greater  volumes  of  air  beinc  punned 
for  the  given  amount  of  heat  transfer  to  take  place. 

Looking  at  the  data  for  the  three  accumulators  shows  that 
the  packed  bed  gives  the  best  compromise  between  heat  exchange 
surface  area  and  volume  of  encapsulation  material.  Another 
factor,  which  must  not  be  overlooked ,  is  flow  distribution. 

In  the  course  of  the  experimental  program  ,  it  was  found  that 
in  the  case  of  the  tube  cross  flow  and  the  tube  parallel  flow, 
all  the  air  flow  tended  to  go  down  one  side  of  the  storage 
unit  due  to  an  elbow  and  a  broad  angled  diffuser  just  upstream 
of  the  storage  unit.  The  packed  bed,  on  the  other  hand  had  the 
advantage  of  evenly  distributing  the  flow  of  air  about  itself. 

V I :  Measurement  System 


The  measurement  system,  consisted  of  two  measurement  chains. 
The  more  elaborate  is  the  thermocouple  temperature  measurement 
chain.  This  starts  with  copper-constantan  thermocouples,  and 
is  shown  in  the  diagram  below: 


The  Mek  II  microprocessor  commands  the  system,  and  the 
multiplexer  scans  U8  channels  of  thermocouple  measurements, 
which  are  then  digitised  by  the  analog-to-digital  convertor 


The  hexadecimal 

n  u  n  b  e  r 

,  ar. 

average  of  1 2 

measurement 

rs  t  V.  o  n 

sequentially  in 

time. 

is  t  r  a 

r.  af'-rr^d  from 

the  memory 

e.  f  t  V.  e 

Mekll  micro r> race 

r»  o r  to  t  h r* 

Pot.  C  o  r.r.  o  d  o  r  ? 

commuter . 

^  V  p  »•  p 

a  software  rrnrrr".  including  the  calibration  curv^  for  th.«» 
t  her  no  c  on  r  1  e  3  r.ibes  the  final  conversion  fro-  the  millivolt 
reading  of  the  'hi  1 1  i  p  1  oxer  to  the  temperature  data  which  are 
required. 

The  second  measurement  chain  is  the  flow  measurement 
chain.  This  is  sim.rly  an  orifice  plate  flow-meter,  with 
a  water  nanometer  used  to  measure  the  pressure  drop  across  it. 
For  this  measurement  ,  a  previously  obtained  calibration  curve 
was  used  to  determine  the  air  mass  flow  rate. 

VII;  Exrie  r  imen  ta  1  Program  and  Results 

The  experimental  program  was  designed  to  test  2  ideas. 

First  of  all,  the  validity  of  the  numerical  model  needed  to 
be  established.  fecond  ,  the  assumptions  going  into  the  model 
needed  to  be  tested.  The  most  dangerous  assumption  is  that 
the  flow  of  heat  into  the  accumulator  is  evenly  distributed 
at  each  point  in  the  cross-section  perpendicular  to  the  flow. 

It  is  a  one-dimer.sional  model.  In  each  of  the  three  accum¬ 
ulators  tested,  the  accumulators  were  instrumented  v'ith  thermo- 
counles  in  the  geometric  center  of  the  cross-section  perpen¬ 
dicular  to  the  flow,  and  at  neveralnoi nts  alonf  the  length  of 
the  accumulators.  This  enabled  a  visualization  of  the  evolution 
of  the  air  and  phase  change  material  temperatures  a  Ions;  the 
length  of  the  bed.  At  three  levels  in  the  accumulators  thermo¬ 
couples  were  put  at  different  points  around  the  cross-section 
in  order  to  test  the  uniformity  of  air  and  phase  change  material 
temperatures . 

A  resume  of  the  experiments  done  follows.  The  primary 
goal  was  to  do  one  complet  charge  and  one  complete  discharge 
of  each  of  the  three  accumulators.  This  was  achieved.  A  large 
amount  of  time  went  into  instrumenting  the  accumulators,  and 
then  finding  the  "bugs"  in  the  system  and  correcting  them. 

There  were  14  tests  altogether,  which  are  tabulated  on 
the  next  page,  together  with  a  list  of  other  oror.aratory  activit 


Date 


Title 


Comments 


Purpose 


2  /  1  / 8 1— 3/14/81  Preparation,  calibration  of  tbernocounler, 


3/iu/ei 

Charge  of 

Tube  Cross 
Flow  unit 

First  charge,  complete  system  vorkec 

3/1  D/61 

ti 

So  me  leakage  of  liquid  PCM  occurei  ; 
flow  distribution  problems 

3/19/81 

Discharge 

0 f  TCF  unit 

Test  terminated  before  complete 
resolidification  had  occured 

3/19/81 

Charge  of 

TCF  unit 

Foam  was  added  for  better  flow 
distribution  in  diffuser 

3/2U/81 

Di schacge 
of  TCF  unit 

8-hour  discharge,  still  not 
sufficient  to  totally  resolidify 

Wo/81 

Charge  of 

TCF  unit 

More  liquid  PCM  leakage  noticed; 
also  leakage  in  ducting  and  dampers 
va s  discovered 

4/10/81 

Discharge  of 
TCF  unit 

Complete  discharge  achieved;  re¬ 
solidification  in  more  than  15  hours 

4/14/81 

Charge  of 

TCF  unit 

t  /  1  4 /8  1 -5  /  1  2 /8 1  a.  Preparation  of  thermocouples 

for  packed  bed  and  tube  parallel  flo- 

b.  Resolution  of  air  leakage  problem 

c.  Mixture  of  PCM  for  packed  bed 
and  tube  parallel  flow  at  Solvay 
chemical  company 

d:  Filling  packed  bed  and  tube 
parallel  flow  units  with  PCM 

5/12/51  Charge  of  Pure  of  leak-proof  air  system; 

TCF  unit  air  flow  measured  downstream  of 

accumulator,  as  well  as  upstream 

5/13/81  Discharge  of  Problems  with  data  acquisition 
TCF  unit  during  test 

5/lW3l  Charge  of  Some  leakage  of  PCM 

Packed  bed 

5/15/31  Discharge  of  No  problems,  complete  discharge 
Packed  bed 

5/19/81  Charge  of  f.'o  PCM  leakage,  very  long  time 

Tube  parallel  required  for  PCM  to  melt 
flow  unit 

7/1/81  Tube  parallel  ?  days  required  for  experiment 
flow  discharge 


Graphs  of  the  experiment al  results  follow:  they  are 
grouped  in  the  three  categories:  tube  cross  flow,  packed 
bed,  and  tube  parallel  flow. 

A.  Tube  Cross  Flow: 

The  tests  on  the  tube  cross  flow  unit  served  to  rive  expe¬ 
rience  with  the  experimental  apparatus,  but  due  to  the  fact 
that  there  was  marked  leakage  in  the  phase  chanre  material, 
they  are  not  conclusive.  Additional  phase  chanre  materia] 
was  not  available,  and  developing  a  system  for  sealing  the 
tubes  at  either  end  would  have  constituted  a  complete  project 
in  itself.  The  problem  lay  in  degradation  of  the  g]ue  at 
either  end  of  the  tubes,  where  plastic  cans  were  fixed. 

In  addition,  since  the  air  flow  rate  going  into  the  accu¬ 
mulator  was  less  than  what  was  measured  at  the  orifice  plate 
upstream,  due  to  air  leaks  in  the  damper  system,  the  convective 
heat  input  was  not  known  with  sufficient  precision  to  warrant 
a  comparison  with  the  numerical  simulation.  A  comparison 
of  the  numerical  results  and  the  experimental  results  for  Test 
number  1  shoved  that  the  phase  chanre  in  the  experiments  proceeded 
much  more  slowly  than  the  computer  simulation  predicted.  This 
was  because  the  mass  flow  rate  which  was  used  as  input  for  the 
simulation  was  that  which  was  measured  at  the  orifice  plate 
upstream  of  the  accumulator,  and  did  not  take  into  account  the 
air  leakage  which  occured  between  the  orifice  plate  and  the 
accumulator  . 

These  limitations  not  withstanding,  some  meaningful 
information  can  be  drawn  from  the  experiments  on  the  tube  cross 
flow  unit.  Graphs  1  and  2  show  the  evolutions  of  the  air  and 
phase  change  material  temperature  s.  These  have  the  characteristic 
plateau  of  fusion.  Also  as  expected,  the  air  temperature  at 
each  point  is  higher  than  the  PCM  temperature  ,  until  steady  state, 
when  the  whole  accumulator  is  at  the  influent  air  temperature, 
is  reached.  The  time  for  reaching  steady  state,  at  60°C  ,  was 
7  hours. 

The  other  information  that  can  he  drawn  regards  the  evenness 
of  the  flow  distribution,  and  the  thermal  performance  as  repeated 
charges  and  discharges  were  performed.  The  second  charre,  as 
shown  in  graphs  3  and  U ,  shows  already  that  the  pretty  plateaus 


1  0 


of  the  first  charge  are  less  flat  than  in  the  first  test. 

Also  the  time  for  chari'inp;  is  about  U  hours,  as  versus  7 
hours  in  Test  1 .  This  should  be  due  both  to  the  fact  that 
scr.e  PCM  had  already  leaked ,  ana  the  fact  that  the  PCM  may  not 
have  been  completely  solidified  when  the  test  was  started. 

The  second  condition  would  mean  that  some  of  the  latent  heat 
of  fusion  remained  in  the  bed  at  the  test's  start. 

Graphs  5,6,  and  7  address  the  nuestion  of  air  temperature 
distribution  about  the  cross  sections  of  the  accumulator. 

These  show  that  in  each  cross  section,  the  temperature  in 
the  center  of  the  tube  on  the  side  (curves  8,11,  and  13) 
were  higher  than  the  temperatures  at  the  2  ends  of  the  central 
tube.  This  can  be  attributed  to  a  higher  heat  capacity  of 
the  greater  amount  of  plastic  at  the  ends  of  the  tubes.  This 
higher  heat  capacity  will  result  in  a  lower  rate  of  temperature 
increase  for  a  given  heat  input.  In  addition  the  air  flow 
rate  should  be  somewhat  higher  at  points  8,11,  and  13  than  at 
points  7 , 9 , 1 0 , 1 2 , 1 U  ,  and  15,  due  to  less  resistance  to  the  flow. 
The  slightly  higher  flow  rate  would  again  result  in  a  higher 
rate  of  temperature  increase  for  the  points  exposed  to  the 
faster  flow. 

Graphs  8,9,  and  10  from  test  number  U  show  that  the  same 
pattern  of  temperature  distributions  about  the  cross-sections 
remains,  even  with  the  additon  of  foam  to  even  out  the  flow 
going  into  the  accumulator.  Thus  one  can  conclude  that  the 
addition  of  the  foam  resulted  in  additional  pressure  loss  which 
was  not  justified  by  an  improvement  in  thermal  performance. 

Tests  number  6  and  7  were  made  to  see  the  temperature  evolution 
during  a  complete  discharge.  Graphs  11,  12,  and  13  show  the 
inlet  and  outlet  temperatures,  their  difference  in  time,  and  the 
integral  of  the  temperature  difference,  respectively.  The 
corresponding  curves  for  the  discharge  are  shown  in  graphs  1 h , 

15,  and  16.  Graph  15  shows  the  phenomenon  of  supercooling 
appearing  in  the  dip  in  the  temperature  difference  and  subsequent 
rise.  Both  graphs  of  the  integrals  of  the  temperature  difference 
should  be  multiplied  by  the  same  mass  flow  rate  to  arrive  at  the 
amount  of  heat  accumulated  or  discharged.  fince  the  interral  is 
not  multiplied  by  mC  .  ir.  either  case, 

*  *  n  i  r* 


the  nur^riri!  valuer 
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cm  be  compared  ;  we  rse  e  that  f  o  r  t  h  °  charge  it  in  1  .11  X  1 r  ^  , 

c 

whereas  for  the  discharge  it  in  3  .  X  1  r,x  .  Thun  on«  third 
of  the  heat  r.ust  be  Lent  through  leakage  .  Irani  17  shown  ve  ]  3 
hov  the  ?(.' "  temperature  follows  the  air  ter.  per  a  t  ur  e  during  the 
course  of  a  charge  of  the  storage  unit,  ar.d  Trap's  1  3  shows  the  same 
s an e  3  curves  for  n  a i s c h ar ce .  graphs  1°,  90,  and  91  show 
that  the  POM  temperatures  are  more  or  less  constant  about,  the 
cross-sections  during  the  discharge. 

The  local  rhe  r.ore  r.on  of  nunercooi  i  nr  can  also  he  seen  fror 
these  graphs.  This  is  a  decrease  of  the  terrnerature  below  the 
fusion  temperature,  followed  by  a  shorn  rise  to  the  fusion  tem¬ 
perature  as  r e so  1 i d i f i c at i o n  occurs,  and  then  a  decrea  e  in  tem¬ 
perature  to  that  of  the  influent  air.  The  supercooling  poses 
problems  for  the  lat ent -heat -of-fus i on  systems  because:  l)  All  the 
heat  is  released  in  a  short  time  when  the  chrystallizatior.  occurs, 
thus  causing  difficulties  in  controlling  the  solidification  and 
2)  During  the  supercooling  ,  air  is  beinr  punned  at  an  expense 
in  electric  energy,  for  no  benefit  in  thermal  energy  transfer. 
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Cranhn  22  through  27  show  the  inlet  and  outlet 
tempera u tre s  ,  their  difference,  and  the  integral  cf  that 
difference  an  a  function  of  tire,  in  the  name  fashion  as  the 
granhs  for  the  charge  and  discharge  for  the  turb  cross  flow. 
Here  the  numerical  valuer,  of  *h-’  ir.terrair,  are  closer  to  each 

r 

other  in  the  charge  and  din  char  re,  at  a.  7  2  X  1  f '  and  h  .  hP 

5 

X  10  ,  respectively. 

Graphs  28  and  20  show  the  air  ar.d  PC”  te^neratur e?  at 
the  top,  middle,  and  bottom  of  th»  accumulators,  and  Grants 
30,31,  and  32  show  the  PCM  temneratures  at  cror.s-'pctions  at 
the  inlet,  middle,  and  outlet  cf  the  accumulator.  It  is  evi¬ 
dent  that  the  temperature  is  very  uniform  across  the  cross- 
sections,  and,  as  mentioned  before,  this  is  a  marked  advan¬ 
tage  of  the  packed  bed.  In  the  discharge  of  the  packed  bed 
there  is  again  the  problem  of  sunercoolinc  ( Graphs  33_3"), 
but  by  having  the  encapsulation  small,  and  thus  nakinr  the 
exit  from  supercooling  occur  in  a  more  dispersed  fashion, 
the  effect  is  smoothed  out  in  the  air  outlet  temerature 
readings.  Also,  the  air  temerature  is  made  more  uniform 
about  the  cross-section  as  the  bed  is  traversed  ( Graphs  36-3*3) 
during  the  charge  of  the  unit. 
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0.  Tube  Pirillel  Flow: 

Hpre  Spain  the  inlet  and  outlet  air  tenoerntures  in 

the  center  of  the  accumulator,  thoir  difference,  and  the 

integral  of  the  difference  are  shown  on  Granhs  3  to  4",  . 

This  can  be  conn a rod  with  confidence  to  the  data  for  the 

racked  bed  charge  ,  since  we  are  confident  of  the  n a  s  r.  flow 

rates.  We  determint  the  average  rate  of  enerpy  storare  as 

follows  :  rt 

Packed  Fed:  rtC  .  (  \AT  -0/t  „ 

n  ,  a  1 r  c  harps 

Jo 

=  .07  k  p  / sec  X  1008  Joule/Ks-°C  X«7?,000  °C-sec 
/ 2 3 0 0 0  sec 
=  1 1f 5 0  .ioule/sec 

Tube  Parallel  flow 

.095  Xp/ sec  X  1008  Joule/Kp-°C 
X  883COO°C-sec  /  70000  sec 
=  1200  joule/sec 

Thus  we  can  see  that  the  packed  bed  stores  heat  at  a  rate 
1.2  times  that  at  which  the  tube  parallel  flow  stores  heat. 

Further  indications  of  the  thermal  performance  of  the 
tube  parallel  flow  unit  are  shown  in  graphs  k'.  and  4«,  .  These 
show  temperature  evolutions  at  different  heiphts  in  the  ac¬ 
cumulator.  Graphs  44  and  45  show  the  inlet  and  outlet  air 
temperature  distributions  about  the  cross-section  of  the 
accumulator.  These  last  2  praphs  show  the  surprisinp 
result  that  the  air  temperatures  are  close  to  each  other  as 
the  accumulator  is  entered,  and  there  is  a  larper  spread  amonr 
them  at  the  exit  from  the  accumulator.  This  may  be  due  to 
irrepular  meltinp  of  the  PCM  about  the  accumulator,  as  well  as 
motion  of  the  melted  PCM  from  one  part  of  the  accumulator  to 
another.  A  complete  discharge  was  performed  on  the  tube  parallel 
flow  unit,  and  the  noteworthy  aspect  of  that  test  war.  that  it 
took  almost  48  hours  for  the  resolidification  to  occur,  because 
of  an  air  inlet  temperature  of  l8°C. 
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VIII:  Concur  i  non  of  numerical  and  Fxperinent  a  1  Results 
Graphs  U6  and  U T  show  the  numerical  curves  and  the 
experimental  points  for  the  racked  bed  unit.  Grach  UC  s  ho  vs 
that  for  the  charge,  the  numerical  simulation  proceeds  faster 
than  the  experimental  data.  A  ft  a  i  n  for  the  discharge,  in  Graph 
i»7,  ti.e  numerical  curves  proceed  more  abruptly  than  experiment. 
This  difference  between  the  numerical  and  experimental  results 
is  probably  due  to  the  resistance  to  heat  transfer  of  the  FCM 
in  the  solid  phase.  Hue h  a  resistance,  neglected  in  the 
numerical  simulation,  would  indeed  result  in  a  slower  rate  of 
discharging  and  char ft  inn.  Thus  if  this  added  resistance  were 
taken  into  account,  the  experimental  and  numerical  results 
would  be  in  better  agreement. 
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I X :  Efficiency  Considerations 

The  mechanical  energy  required  to  store  the  heat  in  each 
accumulator  can  be  vritten  as  4!’  X  V  =  Power  =  V .  Thin  nays 
that  the  power  required  to  maintain  a  volume  flow  of  v  across  a 
pressure  difference  of  AP  ,  is  equal  to  the  product  of  the  two. 
Looking  at  this  in  terms  of  work,  it  car.  be  said  that  W  is  the 
amount  of  work  done  in  moving  a  volume  of  fluid  V  across  a 
pressure  drop  A?.  Thus  we  have  the  amount  of  work  required  for 
our  accumulators  by  measuring  the  pressure  drop  ac cross  them. 

The  thermal  energy  transfer  is  mC  .  AT,  vhereAT  is  the 

p,  air 

difference  between  the  inlet  and  outlet  air  temperatures,  and 
m  is  the  mass  flow  rate.  Ve  then  find  Q  =  pvC  A  T.  Now  if 
we  want  a  measure  of  the  amount  of  heat  transferred  in  relation 
to  the  amount  of  work  required  to  transfer  that  heat,  we  can 
devide  Q  by  W.  This  is  simply  C^AT/AP.  This  then  is  readily 
available  data  from  the  experiments,  and  it  has  been  found  that 
in  every  case  it  is  on  the  order  of  50  to  100.  This  high 
apparent  efficiency  though  must  be  considered  to  be  divided 
first  by  2,  since  there  is  a  charge  and  a  discharge  of  the 
accumulator  in  any  application  where  the  heat  is  used.  Then  it 
mus‘  again  be  divided  by  3,  for  the  overr.ll  efficiency  of 

the  electric  system.  This  brings  it  down  bv  almost  an  order 
of  magnitude,  to  8  to  l6  already.,  down  from  50  to  100. 

Together  with  this  idea  of  efficiency,  a  criterion  for 
stopping  the  charge  or  discharge  should  be  developed.  That  is, 
the  W  remains  more  or  less  constant  during  the  charge  and  dis¬ 
charge  while  Q  decreases  with  time  due  to  AT  decreasing  with 
time.  Thus  the  efficiency  is  decreasing  throughout  the  experi¬ 
ment  and  after  a  certain  point,  oil  or  electricity  could  be  used 
more  advantageously  than  by  pumping  low  temperature  air  in  and 
out  of  a  storage  unit. 


i 
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x :  coi!CLU3io::s  a::d  rkcommohdatioiis 

The  experimental  study  shows  that  from  a  heat  transfer 
standpoint  the  packed  bed  in  the  best  method  of  encapsulation. 
It  has  a  negligible  pressure  loss,  large  heat  exchange  sur¬ 
face  area,  a  self  distribution  of  the  flow  and  an  avoidance 
of  the  problems  of  supercooling.  On  the  other  hand,  the  geo¬ 
metry  which  was  worst  in  heat  transfer  characteristics ,  the 
tube  parallel  flow,  was  the  only  one  of  the  three  that  had  no 
PCM  leakage  problem,  while  the  packed  bed, though  the  container 
were  sealed  with  silicon  rubber  on  the  caps,  did  have  leakage 
problems.  Thus  the  problem  is  one  of  feasibility  and  cost  of 
construction  as  well  as  heat  transfer. 

Further  studies  of  such  storage  units  should  include  an 
operation  over  a  broau  range  of  Reynolds  numbers,  to  find  the 
optimum  range  for  thermo-mechanical  efficiency.  Also, 
experiments  should  be  done  which  include  solar  collector  w: 
inlet  temperature  varying  with  time,  instead  of  the  fixed 
inlet  temperature  which  was  used  in  the  past  experiments. 
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Photographs 


Appendix  I 

of  Experimental  Program 
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